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ABSTRACT 

We present high-resolution images of the Sun in the near ultraviolet spectral range between 214 nm 
and 397 nm as obtained from the first science flight of the 1-m Sunrise balloon-borne solar telescope. 
The quiet-Sun rms intensity contrasts found in this wavelength range are among the highest values ever 
obtained for quiet-Sun solar surface structures - up to 32.8 % at a wavelength of 214 nm. We compare 
with theoretical intensity contrasts obtained from numerical magneto-hydrodynamic simulations. For 
388 nm and 312 nm the observations agree well with the numerical simulations whereas at shorter 
wavelengths discrepancies between observed and simulated contrasts remain. 
Subject headings: Sun: general — Sun: atmosphere — Sun: granulation 



1. INTRODUCTION 

Quiet-Sun intensity fluctuations provide fundamental 
information on the thermal structure of the convective 
overshoot region at the solar surface. A large number 
of measurements of the quiet-Sun intensity fluctuations 
have been made in the visible spectral ra n ge, st arting 
in the 1950s (jFrenkiel & Schwarzschildl 1 19551) and 
reaching to comparisons between recent seeing-frce 
observations from the Solar Optical Telescope (SOT) 
on the Hinode satellite, which provide contrast values 
comparable to t hose from state - of-the - art numerical 
simul ations fe.g.. | Danilovic et al.l [20 08: Mathe w et al.| 
2009; iWedemever-Bohm fc Rouppe van der Voortl 

2009). Overviews of the results of former me asure- 
ments are given by Beckers fc ParnelH (|1969f ) and 
ISanchez Cuberes et al.1 (|2000h . Although many intensity 
contrast measurements are based on data obtained after 
careful correction of atmospheric and/or instrumental 
effects, it turns out that insufficient knowledge of optical 
aberrations and stray light still are the major obstacles 
for accurate contrast measurements. 

In the near ultraviolet (NUV), the intensity fluctua- 
tions of the solar surface region have been hitherto largely 
unknown, owing to the lack of reliable imaging data. The 
strong dependence of the Planck function in the NUV to 
the temperature makes the intensity highly sensitive to 
temperature fluctuations. Therefore, intensity contrasts 
also represent important diagnostics to probe the validity 
of numerical simulations of convection at granular scales. 

From an observational point of view, it is well known 
that seeing effects and straylight contributions increase 
towards shorter wavelengths. Ground-based NUV ob- 
servations are, additionally, hampered by strong at- 



mospheric absorption (mainly by stratospheric ozone). 
Below the atm ospheric cut-off at about 315 nm (see 
Gandorfer 2005) , observations cannot be carried out from 
the ground. 

In the present study we show for the first time high- 
resolution imaging data and intensity fluctuations in 
the NUV region down to 214 nm. These data were 
obtained during the first science flight o f the balloon- 
borne 1-m SUNRISE observato r y (see iBarthol et al 
20101: iMartmez Pillet et~aT1 120101: fBerkefeld et aT ~ 
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Solanki et al.l l2010f T We present disk-center quiet-Sun 



hirzbcrger@mps . mpg.de 

* The National Center for Atmospheric Research is sponsored 
by the National Science Foundation 



NUV intensity contrast s obtained with the Su nrise Fil- 
ter Imager (SuFI, see iGandorfer et al.1 120101 ) onboard 
Sunrise. 

2. observations 

The first science flight of the balloon-borne Sunrise 
telescope was carried out between June 8 and June 13, 
2009, from Kiruna, Sweden, to Somerset Island, Canada. 
The average cruise altitude was about 36km, i.e., at 
the upper boundary of the stratospheric ozone layer, 
so that high resolution images in the NUV could be 
obtained. The spectral regions observed by SuFI are 
214 nm, 300 nm, 312 nm (close to the bandhead of the 
OH molecule), 388 nm (bandhead of the CN molecule), 
and 396.8 nm (CallF line core). The widths (FWHM) 
of the corresponding SuFI filters are 10 nm, 5nm, 1.2 nm, 
0.8 nm, and 0.18 nm, respectively. The passbands within 
CN and OH bandheads were selected because of their 
high temperature sensitivity, especially UV OH lines 
which are favorable f or imaging of both the quiet photo- 
sphere and sunspots ([Berdyugina et aLll2003[ ). The filters 
are mounted in a filter wheel that allows for sequential 
observations of the five wavelength channels. 

SuFI is equipped with a 2k x 2k CCD detector (Pix- 
elVision BioXight) at an effective focal length of 121 m. 
The camera pixel pitch is 12 /xm corresponding to 0"0207 
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(at 300 nm) on the sky. For post-facto correction of 
low-order optical aberrations by means of phase-diversity 
wavefront sensing (PD), SuFI was equipped with a de- 
vice providing a focussed and a defocussed image, each 
on one half of the detector. The effective field of view 
(FOV) is thus limited to about 15" x 40". 

The exposure times for the SuFI observations were 
strongly wavelength-dependent and ranged from 75 ms 
in the CN channel up to 30 s in the 214 nm channel. In 
this latter channel, the intensity level was high enough 
for detection only around maximum solar elevation and 
maximum balloon altitude. Therefore, observations in 
the 214 nm channel were only carried out around local 
noon. 

Here we show data obtained on June 11, 2009, be- 
tween 20:00 UT and 21:00 UT (300 nm, 312 nm, 388 nm, 
and 397 nm) and on June 9, 2009, between 14:00 UT and 
15:00 UT (all SuFI channels). Examples of single images 
obtained in the 5 SuFI wavelength channels are shown in 
Fig.ffl 

3. DATA ANALYSIS 

The raw ("level 0") SuFI images were corrected for 
dark and flat fields. Dedicated flat-field images were ob- 
tained by moving the telescope during exposure. For flat 
field correction of the science images accumulations of 
all flat-field and science images obtained in an interval of 
±2 hours around the observation time were used. In ad- 
dition, clusters of bad detector pixels were eliminated by 
applying a local median filter. Residual inhomogencitics 
originating from scratches at the entrance window of the 
SuFI camera were removed by local smoothing using low- 
pass filters. 

The flat- and dark-field corrected SuFI images (hence- 
forth labelled as "level 1" data) were reconstructed by 
means of a PD al gorithm. We used a PD c ode origi- 
nally described by lLofdahl fc Scharmerl (11994ft and fur- 
ther developed as d escrib ed by Bonet et al.l (|2004l) and 
iVargas Dommgue3 (|2008f l. In order to optimize the re- 
construction of the SuFI data, several additional features 
were included, which will be described together with a 
performance study of the new PD c ode in a forthcoming 
publication (|Hirzberger et al.ll2010L in prep.). 

In contrast to ground-based observations, the contri- 
bution of atmospheric influences to the total wavefront 
aberrations at balloon altitudes is small compared to 
instrumental effects. Therefore, we expect isoplanatic 
conditions, i.e. constant wavefront aberration across the 
SuFI FOV. Temporal variations of the wavefront defor- 
mations are assumed to arise only from changes of the 
instrument's temperatures and from slowly varying me- 
chanical deformations, which both are mainly correlated 
with the solar elevation. These assumptions allow a ho- 
mogeneous reconstruction of the SuFI images. We thus 
included a procedure, which permits averaging the fit- 
ted Zernike coefficients across the entire FOV as well 
as over all images within a certain time span (usually 
one hour). Subsequently, all images taken during this 
time span were reconstructed with constant and highly 
reliable wavefront errors. PD reconstruction of individ- 
ual images is, henceforth, denoted as "level 2" reduction, 
while images reconstructed with averaged wavefront er- 
rors are labelled as "level 3" . 



4. RESULTS 

4.1. Disk-center intensity contrasts 

The quiet-Sun disk-center rms intensity contrasts, 
filrms, derived from the data obtained on June 11, 2009 
are plotted in Fig. [2l In each spectral channel, for 
180 consecutive images (approximately 24min) the im- 
age stabilization syst em (Sunrise Correl ating Wave- 
front Sensor, CWS; see lBerkefeld~et~a!1 l2010') was locking 
in closed loop. The obtained intensity contrasts show 
considerable temporal fluctuations which are partially 
caused by the evolution of the solar structures and p- 
modes but also by a variation of the image quality. The 
image quality is affected, first, by the limited overall 
stability given by the gondola pointing and the image 
stabilization in presen ce of wind gusts and vibrations 
(|Berkefeld et al.l [2010D . and, secondly, by residual arti- 
facts, such as scratches and bad pixels in the level- 1 im- 
ages (see Sect. [3]), which are amplified by the PD recon- 
struction. After visual inspection, the influence of the 
latter is considered to be negligible in the data presented 
here. 

In addition to the temporal fluctuations, a rather con- 
stant offset of approximately 2 — 3 % (depending on wave- 
length) of the intensity contrasts from level-2 to level-3 
data is evident. This offset is expected since for the 
level-2 data always the local best fits to the wavefront 
aberrations are used for reconstruction (which also may 
cause local over-reconstructions), whereas the averaged 
wavefronts used for reconstructing the level-3 data may 
locally underestimate the wavefront errors (e.g. due to a 
small amount of anisoplanatism). 

Irrespective of the effect of the image processing meth- 
ods, the photospheric intensity contrasts shown in Fig. [2] 
are among the highest quiet-Sun contrasts ever mea- 
sured. The contrast variations within the time series 
are to a lesser extent related to granular evolution than 
to the varying image quality due to residual pointing 
errors. This sometimes leads to artifacts due to local 
over-reconstruction and thus to unrealistically high con- 
trast values. In the majority of images this effect, how- 
ever, tends to result in reduced intensity contrasts (by 
smearing out solar surface structures). In order to al- 
low for this systematic reduction of contrasts and to clip 
outliers, we consider the "mean-plus-one-sigma" values, 
{5I rms )+(T, as reliable maximum contrast measurements. 
Here, (5I rms ) denotes the temporally averaged rms con- 
trast and a denotes the corresponding standard devia- 
tion. The mean-plus-one-sigma values are overplottcd in 
Fig. [2] and given in Tab. Q] 

Figure [3] shows the intensity contrasts of 15 images 
from the 214 nm and 300 nm channels of the June 9 data. 
As a result of the long exposure time at 214 nm, the 
number of reliable images obtained during the observing 
period is limited to that small number. The mean-plus- 
one-sigma rms contrasts of the June 11 data are marked 
by arrows. Due to higher residual pointing errors, which 
caused some smearing of the images, the June 9 data are 
of slightly worse quality than the June 11 data. This 
particularly affects the 214 nm data. Nevertheless, the 
obtained intensity contrasts in the 214 nm channel are 
the highest photospheric quiet-Sun contrasts ever mea- 
sured. 
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4.2. Intensity contrasts from simulations 

In order to compare with the measured intensity con- 
trasts, we have calculated NUV intensity maps from nu- 
merical magneto-hydrodyn amic (MHD) simu lations of 
the quiet solar photosphere (Vo gler et al.ll2005l ). We used 
snapshots from a simulation run with a mean vertical 
magnetic field strength of (B z ) = 50 G and a horizon- 
tal cell size of 20.8 km. The intensity maps were cal- 
culated by us ing the spectral synthesis code ATLAS9 
(|Kuruczl fl993) whi ch specifies the opacit y distribution 
functions (ODFs, iStrom fc Kuruczl 1 19661 ) in the SuFI 
300 nm, 312 nm, and 388 nm channels. With the cur- 
rent implementation for computing ODFs, we are un- 
able to calculate reliable intensities below a wavelength 
of 220 nm. In order to estimate the simulated contrast at 
214 nm, we tentatively applied the ODF for 220 nm to the 
MHD snapshots. In addition, in the SuFI 388 nm channel 
a full spectral synthesis wi t h the SPINOR code rtSqlankil 
Il987t iFrutiger et all l2000t iBerdvugina eTaTl 12003ft .~h> 
eluding all available atomic and molecular line parame- 
ters within the SuFI 388 nm spectral bandpass have been 
performed. 

Nine synthetic intensity images at 300 nm and the 
corresponding SI rrns for all SuFI channels are shown 
in Fig. |4j Owing to the small area (6 Mm x 6 Mm) of 
the MHD box, the simulated contrasts show a temporal 
variation with a standard deviation of about 1 % (ex- 
cept for the 220 nm spectral region where the value is 
about 3 %) . The contrasts obtained with SPINOR in the 
388 nm channel are consistently about 1.3% higher than 
those obtained by using ODFs. This small difference, 
comparable to the temporal variation of <5/ rms , indicates 
that the ODF-based spectral synthesis in the other spec- 
tral channels may not be completely unrealistic, although 
NLTE effects have not been taken into account. The 
mean contrasts obtained from the MHD simulations are 
also given in Tab. [TJ 

5. DISCUSSION AND CONCLUSIONS 

The NUV intensity contrasts obtained from Sun- 
RISe/SuFI data are systematically lower than the values 
from MHD simulations by several percent (but nearly a 
factor of two at 214 nm). Since the spatial resolution 
of both datasets is similar and because the instrumental 
aberrations have been removed by means of PD recon- 
struction, three possible sources for the discrepancy re- 
main: (i) inaccurate physics in the simulations, (ii) omis- 
sion of NLTE effects in the intensity calculations, (iii) 
scattered light. Before we can estimate the importance 
of NLTE effects or the need to improve the simulations, 
we need to judge the influence of scattered light on the 
measured contrasts. 

The amount of scattered light in the SuFI data has not 
yet been fully assessed. A first estimate for the 300 nm, 
312 nm and 388 nm SuFI channels was derived from limb 
observations. The straylight components of the modula- 
tion transfer function obtained from limb profiles at the 
different wavelengths were applied to the numerical sim- 
ulations to estimate the effect of straylight on the rms 
int ensity contr a sts (s ee Tab. Q] details will be presented 
bv lFeller et all (|20ToL in prep.)). For 214nm and 397nm 
we were not able, yet, to obtain a straylight estimate due 
to the lack of limb observations (214 nm) and due to the 



ubiquitous presence of spicules (397 nm). At 312 nm and 
388 nm, the synthetic images after stray light contamina- 
tion display contrast values that lie between the level- 2 
and the level-3 data. At these wavelengths, the simula- 
tions thus give results consistent with the measurements. 
At 300 nm the simulations give a value of 25.5% that is 
still 1.3 % higher than the 24.2 % obtained from the level- 
2 data. 

Estimates of scattered light contributions in 
(also seeing free) data obtained with Hinod e /SOT 
were carried out by IMathew et al.l (120091) an d 
IWedemever-B5hm fc Rouppe van der Voortl ((2009). 
They found a significant increase of the intensity 
contrasts in the visible Hinode channels, after decon- 
volving with point-spread functio ns for scattere d light . 
For the 388 nm spectral region IMathew et all {2009) 
found a quiet-Sun intensity contrast of 21.8% after 
removing scattered light, which is very similar to the 
mean-plus-one-sigma value of our level-2 data from 
June 11, 2009 (cf. Tab. [J). Therefore, we conclude that 
after applying straylight correction to the SuFI-388nm 
da ta, the obtaine d cont rasts will exceed those obtained 
by IMathew etail (|2009f ) fr om Hinode data and will be 
comparable to synthetic contrasts achieved from MHD 
data. 

The consistency between the measured and the nu- 
merically simulated contrasts indicates that the tem- 
perature fluctuations in the lower photosphere are cor- 
rectly described by the hydrodynamical simulations, 
which predict values of 5T rms between 2 % and 5 % 
on surfaces of constant optical depth (at 500 nm) be- 
tween T500 = 1 and T500 = 0.01. Taken together with 
the repro duction of shifts and asymmetries of spectral 
lines Ce.g.lNordlund etall l2009) and inverse granulation 
([Cheung et al.l 120071 ) bv simulation results, our contrast 
data thus provide evidence for an extension of the con- 
vective overshoot to a height of roughly 300 km above 
the average level of optical depth unity, at which height 
the rms fluctuations of the vertical velocity reach a lo- 
cal minimum after a steep decrease (by about a factor of 
4) from their maximum slightly below the visible surface 
(optical depth unity). 

Direct translation of rms contrasts into e.g. tempera- 
ture stratification of the solar photosphere can be car- 
ried out only by means of simplified models estimat- 
ing the brightness tempera t ure of atmospheric features 
(e.g. lSolanki fc Unruhlll99l iSobotka et al.ll2000D . In the 
NUV the rms intensity contrast of broad band images as 
obtained with Sunrise/SuFI is, however, only partly a 
function of height. It depends just as much on the tem- 
perature sensitivity of the Planck function and of the 
lines in the SuFI passbands, the atomic species, ioniza- 
tion stages and in the case of molecular transitions, the 
dissociation energy of the molecules. 

In summary, we conclude that the data obtained from 
the first science flight of Sunrise/SuFI are of resounding 
quality. They show the highest quiet-Sun photospheric 
intensity contrasts ever measured, even without correct- 
ing for scattered light. 
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TABLE 1 

MEAN-PLUS-ONE-SIGMA VALUES OF SIrms (SEE FlGS. [2] AND [3j) AND MEAN 5I r ms RESULTING FROM MHD SIMULATION DATA (IN PERCENT). 
IN THE RIGHTMOST TWO COLUMNS, MEAN SIrms OF MHD DATA CONVOLVED WITH PRELIMINARY ESTIMATED LEVELS OF SCATTERED LIGHT 

ARE GIVEN. 
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Fig. 1. — Examples of PD-reconstructed (level 2) images obtained in the 5 Sunrise/SuFI wavelength channels. The displayed image in 
the leftmost panel (214 nm) was obtained on June 9, 2009, 14:59:24 UT, the other images on June 11, 20:07:41 UT (300 nm), 20:07:42 UT 
(312nm), 20:07:38 UT (388 nm), and 20:07:39 UT (397 nm). All images show fields of quiet granulation near disk center. The intensities 
have been normalized to the mean intensity I over the FOV and the scales have been limited to / ± 10cr (214 nm) and J ± 5<r (all other 
channels), respectively, where a denotes the standard deviation of the corresponding intensity. 
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Fig. 2. — rms intensity contrasts of 180 images in each of four wavelength channels, obtained on June 11, 2009, between 20:00 and 
21:00 UT. Solid lines correspond to level-2 data and dotted lines refer to level-3 data. The dashed and dash-dotted horizontal lines lie one 
standard deviation above the mean contrast of the level-2 and level-3 images, respectively. In the 397 nm band, the contrast is slightly 
higher in the second half of the image series than in the first half since the FOV had moved by a few arcseconds on the solar disk. 
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Fig. 3. — rms intensity contrasts of 15 images taken in the 214 nm (panel a) and 300 nm (panel b) wavelength channels, obtained on June 9, 
2009, between 14:00 and 15:00 UT. Solid lines correspond to level-2 data and dotted lines refer to level-3 data. The dashed and dash-dotted 
horizontal lines show the mcan-plus-onc-sigma values of level-2 and level-3 images, respectively. For comparison, the mean-plus-one-sigma 
values obtained from the June 11 300 nm data are indicated by arrows. 
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Fig. 4. — Upper panel: Nine brightness images at 300 nm from a 
MHD simulation with a mean vertical field strength of (B z ) = 50 G; 
lower panel: corresponding <5/ rms calculated for four of the five 
SuFI spectral channels, using two different methods for 388 nm 
(ODF- based and full spectral synthesis). Since ODFs do not allow 
reliable intensity calculations at 214 nm, the computations have 
been carried out at the nearest wavelength at which they are better 
established, 220 nm. 



